THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING 



PLAMEN STEFANOV AND GUNTHER UHLMANN 

Abstract. We study the mathematical model of thermoacoustic and photoacoustic tomography 
when the sound speed has a jump across a smooth surface. This models the change of the sound 
speed in the skull when trying to image the human brain. We derive an explicit inversion formula 
in the form of a convergent Neumann series under the assumptions that all singularities from the 
support of the source reach the boundary. 



1. Introduction 

In this paper, we study the mathematical model of thermoacoustic and photoacoustic tomography 
for sound speed that jumps across a smooth closed surface. The physical setup is the following. A 
short impulse of microwaves or light is sent through a patient's body. The cells react by emitting 
an acoustic signal that is being detected on a smooth surface around the patient's body. Then the 
problem is to recover the density of the source of the acoustic waves that can be used to recover 
the absorption rate of the tissue at any point, thus creating an image, see e.g., [H |9j [HI [T5l [33J, 
126] . For more detail, an extensive list of references, and the recent progress in the mathematical 
understanding of this problem in the case of a constant or a smooth sound speed, we refer to 

m mini nana m- 

Let Q C R n be a bounded domain with smooth boundary. Let r C Si be a smooth closed, 
orientable, not necessarily connected surface. . Let the sound speed c(x) > be smooth up to V 
with a nonzero jump across it. For x G T, and a fixed orientation of T, we introduce the notation 

(1-1) c int (x) = cL , c cxt (x) = cL 

1 -L mt 1 L ext 

for the limits from the "interior" and from the "exterior" of O \ T. Our assumption then is that 
those limits are positive as well, and 

(1.2) c int (x) / Cext(aO, Vx G T. 

This problem was proposed by Lihong Wang at the meeting in Banff on inverse transport and 
tomography in May, 2010 and it arises in brain imaging [34, 35J. In that case, the brain is represented 
by some domain f^o d f2. Let f2i be another domain representing the brain and the skull, so that 
<s ili (s and Oi \ On is the skull, see Figure [l] The measuring devices are then typically 
placed on a surface encompassing the skull, modeled by dQ in our case. Then 

with the speed jumping by about a factor of two inside the skull fii \ Qq. Another motivation 
to study this problem is to model the classical case of a smooth speed in the patient's body but 
account for a possible jump of the speed when the acoustic waves leave the body and enter the 
liquid surrounding it. 
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The mathematical model can then be described as follows. Let u solve the problem 

in (0,T) x R n , 



' (9 t 2 -c 2 A)u 


= 


u \r- 


1 -L CXt 


du 1 


du 1 




~~ 9" 1 Text 


u\t=o 


= /, 


d t u\ t =o 


= o, 



(1.3) 



where T > is fixed, u|r intcxt is the limit value (the trace) of u on T when taking the limit from 
the "exterior" and from the "interior" of V, respectively. We similarly define the interior /exterior 
normal derivatives, and v is the exterior unit (in the Euclidean metric) normal to T. 

Assume that / is supported in Cl, where f2 C R n is some smooth bounded domain. The mea- 
surements are modeled by the operator 

(1-4) Ai/ := u\ [0jT]xdn . 

The problem is to reconstruct the unknown /. 

We denote more generally by A[/i, fy] = Ai/i + A2/2 the measurements corresponding to general 



Cauchy data [/i,/2] in (1.3). We will work with fx, f2, supported in some compact TC in f2. In 
applications, this corresponds to f\, /2, that are not necessarily zero outside 1C but are known there. 
By subtracting the known part, we arrive at the formulation that we described above. We also 
assume that c = 1 on R n \ $7. We formulate the main results for Cauchy data [/i,0] for simplicity 
of the exposition but we do most of the preparatory work for general Cauchy data. 

The propagation of singularities for the transmission problem is well understood, at least away 
from possible gliding rays [10} 129] 120] [21] . When a singularity traveling along a geodesic hits the 
interface V transversely, there is a reflected ray carrying a singularity, that reflects at V according to 
the usual reflection laws. If the speed on the other side is smaller, there is a transmitted (refracted) 



ray, as well, at an angle satisfying Snell's law, see (4.43). In the opposite case, such a ray exists 
only if the angle with V is above some critical one, see (4.44). If that angle is smaller than the 
critical one, there is no transmitted singularity on the other side of V. This is known as a full 
internal reflection. This is what happens in the case of the skull when a ray hits the skull boundary 
from inside at a small enough angle, see Figure [TJ Therefore, the initial ray splits into two parts, 
or does not split; or hits the boundary exactly with an angle equal to the critical one. The latter 
case is more delicate, and we refer to section [4] for some discussion on that. Next, consider the 
propagation of each branch, if more than one. Each branch may split into two, etc. In the skull 
example, a ray coming from the interior of the skull hitting the boundary goes to a region with 
a smaller speed; and therefore there is always a transmitted ray, together with the reflected one. 
Then a single singularity starting at time t = until time t = T in general propagates along a few 
branches that look like a directed graph. This is true at least under the assumption than none of 
those branches, including possible transmitted ones, is tangent to the boundary. 

If /2 = 0, that is the case we are interested in, singularities from (xo,£°) start to propagate in 
the direction £° and in the negative one if none of the branches reaches d£l at time T or 

less, a stable recovery is not possible [23]. In section [2j we study the case where the initial data 
is supported in some compact K. C ft \ T and for each (xo,£°) £ T*/C \ 0, each ray through it, 
and through (a?o,£°) has a branch that reaches 50 transversely at time less than T. The main 
idea of the proof is to estimate the energy that each branch carries at high energies. If there is 
branching into non-tangent to the boundary rays, we show that a positive portion of the energy 
is transmitted, and a positive one is reflected, at high energies. As long as one of these branches 
reaches the boundary transversely, at a time at which measurements are still done, we can detect 
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Figure 1 . Propagation of singularities for the transmission problem in the "skull" 
example. The shaded region represents the "skull", and the speed there is higher 
than in the non-shaded part. The dotted curves represent the propagation of the 
same singularity but moving with the negative wave speed. 



that singularity. If we can do that for all singularities originating from /C, we have stability. This 
explains condition (2.9) below. Uniqueness follows from unique continuation results. 

Similarly to [26], assuming (2.9), we also get an explicit converging Neumann series formula for 

As in the case of a smooth speed considered in (26j the "error" 



assuming (2.9 
reconstructing f , see Theorem 



2.1 



operator K in (2.10) is a contraction. An essential difference in this case is that K is not necessarily 
compact. Roughly speaking, Kf corresponds to that part of the high frequency energy that is still 
held in Q until time T due to reflected or transmitted signals that have not reached dQ yet. While 



the first term only in (2.10) will still recover all singularities of f, it will not recover their strength, 
in contrast to the situation in [22], where the speed is smooth. Thus one can expect somewhat 
slower convergence in this case. 

Acknowledgment. We would like to thank Lihong Wang for helpful conversations on brain 
imaging using thermoacoustic tomography. 



2. Main result 



Let u solve the problem 



(2.1) 





= 




= it L 


au 


du 1 




~~ du\ 


u\ t =Q 


= h, 


< dtu\ t =o 


= f2, 



in (0, T) x R" 



where T > is fixed. Compared to (1.3), we allowed /2 to be non-zero but in the main result, we 
will take /2 = for simplicity of the exposition. Set f = [/i, /g]. Assume that / is supported in Vt, 
where f2 C R n is some smooth bounded domain. Set 



(2.2) 



Af 



u \[0,T]xdQ- 
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The trace Af is well defined in C(o)([0, T]; H l l 2 (dVL)}, where the subscript (0) indicates that we 
take the subspace of functions h so that h = for t = 0. For a discussion of other mapping 
properties, we refer to [13], when c has no jumps. By finite speed of propagation, one can reduce 
the analysis of the mapping properties of A to that case. 
In the thermoacoustic model, fi = 0. For this reason we set 

(2.3) Ai/i:=A[/ x ,0]. 

This notation is justified by setting Ai 2 to be the components of A (that sends vector functions to 
scalar functions, i.e., Af = A1/1 + A2/2. In this paper, we use boldface to denote vector functions 
or operators that map scalar or vector functions to vector functions. 

The standard back-projection that would serve as some kind of approximation of the actual 
solution is the following. We cut off smoothly Ai fx near t = T to satisfy the compatibility conditions 
in the next step; and then we solve a backward mixed problem with boundary data the so cut Ax fx) 
and Cauchy data [0, 0] at t = T. As in the case of a smooth speed, see [114 126], one can show that 
such a back-projection would converge to /, as T — > 00 at a rate that depends of /; and at least 
at a slow logarithmic one, if one knows a priori that / G H 2 , see pi . If T = dti , where O C O 
is strictly convex, then in the case that the speed outside fio is faster than the speed inside (then 
there is full internal reflection), the convergence would be no faster than logarithmic, as suggested 
by the result in |22j . In the opposite case, it is exponential if n is odd, and polynomial when n is 
even [3J. Our goal in this work is to fix T however. 

In [SB], we proposed the following modified back-projection. Given h that we will be chosen to 
be Ax fx later, let v solve 

' (d t 2 - c 2 A)v = in (0, T) x ft, 
(2 4) J v \[o,T]xdn = K 

I v\ t =T = <P, 

d t v\t=T = 0, 

where <j) solves the elliptic boundary value problem 

(2.5) A0 = O, <f)\ an = h(T,-). 
Then we define the following pseudo-inverse 

(2.6) Ah:=[v(0,-),v t (0,-)]=:[Axh,A 2 h] in 0. 
By [IT] , and using finite speed of propagation [32] , one can show that 

A : Hf Q) {[0,T\ x dfl) -> U ^ H%(Sl) x L 2 (n) 

is a continuous map. Note that the mapping properties above allow us to apply A to Af only when 
f is compactly supported in Q but the theorem above shows that AA extends continuously to the 
whole The function Axh with h = Ax fx can be thought of as the first "approximation" of 

fx- On the other hand, the proof Theorem |2.1| below shows that it is not a good approximation, 
see Remark 12. II 

To explain the idea behind this approach, let us assume for a moment that we knew the Cauchy 
data [it, ut] on {T} x S7. Then one could simply solve the mixed problem in [0, T] X with that 
Cauchy data and boundary data A\f\. Then that solution at t = recovers /. We do not know 
the Cauchy data [it, ut] on {T} x O, of course, but we know the trace of u (a priori in H 1 for t 
fixed) on {T} x <9f2. The trace of ut does not make sense because the latter is only in I? for t = T. 
The choice of the Cauchy data in (|2.4[) can then be explained by the following. Among all possible 
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Cauchy data that belong to the "shifted linear space" (the linear space H(£l) translated by a single 
element of the set below) 

(2.7) {g=[g ll g 2 }€H 1 (n)(BL 2 (n); gi \ dn = h(T, •)} 

we chose the one that minimizes the energy. The "error" then will be minimized. We refer to 
Lemma 5.1 and the proof of Theorem 2.1 for more details. 

The discussion in the Introduction and in section [4] indicates that the singularities that we are 
certain to detect at d£l lie in the following "non-trapped" set 

U =Ux,£) G S*(Q \ T); there is a path of the "geodesic" issued from either 

(2.8) 1 V ' - 

(x, £) or (x, — £) at t = never tangent to T, that is outside at time t = Tj. 

Actually, U is the maximal open set with the property that a singularity in IA "is visible" at 
[0, T] x d£l; and what happens at the boundary of that set, that includes for example rays tangent 
to r, will not be important for our analysis. We emphasize here that "visible" means that some 
positive fraction of the energy and high frequencies can be detected as a singularity of the data; 
and of course there is a fraction that is reflected; then some trace of it may appear later on <9f2, 
etc. 

One special case is the following. Take a compact set K, C f2 \ T, and assume that 

(2.9) S*KdU. 

In other words, we require that for any x G K, and any unit £ G S*tC, at least one of the multi- 
branched "geodesies" starting from (x,^), and from (x, — £), at t = has a path that hits d£l for 
time t < T and satisfies the non-tangency assumption of (2.8). Such a set may not even exist for 
some speeds c. 

Example 1. Let flo C be two concentrating balls, and let c be piecewise constant; more 
precisely, assume 

n = B(0,R), n = B(0,R ), 0<R <R, 

and let 

c = Co < 1 in f^oi c = 1 in R n \ Qq. 

Then such a set fC always exist and can be taken to be a ball with the same center and small enough 
radius. Indeed, the requirement then is that all rays starting from fC hit T at an angle greater than 
a critical one n/2 — ocq, see (4.44). This can be achieved by choosing fC = B(0, p) with p Rq. An 
elementary calculation shows that we need to satisfy the inequality p/Ro < sinao = Co, i.e., it is 
enough to choose p < cqRo- Then there exists To that is easy to compute so that for T > To, (2.9) 
holds. If c = Co > 1 in f^o, then any compact K, in satisfies (2.9). In that case, there is always a 
transmitted ray leaving S7o- 

Example 2. This is a simplified version of the skull model. Let fio C Oi C fi be balls so that 

Q = B(0,R), n = B(0,R ), Oi = B(0,i?i) 0<R <R 1 <R, 

Assume that 

c|n = c 0, c lfii\fi = c l> c Ir"\^i = 1 
with some constants Co, c\ so that co < c±, c\ > 1. Here, Co models the acoustic speed in the 
brain, c\ is the speed in the skull, and 1 is the acoustic speed in the liquid outside the head. If 
for a moment we consider VLq and only, we have the configuration of the previous example. If 
K. = -6(0, p) with p < (co/ci)i?o, then K, satisfies (2.9). Now, since c\ > 1, rays that hit dfti always 
have a transmitted part outside and therefore (2.9) is still satisfied in fi. 
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Let ILjc : Hd(£1) — > Hd(1C) be the orthogonal projection of elements of the former space to the 
latter (considered as a subspace of Hd(£1)). It is easy to check that Hfcf = f\x. — PdK(f\dK.)i 
where Pqjq is the Poisson operator of harmonic extension in /C, see also Lemma |5.1[ Let ljc be the 
restriction of functions defined on O, to /C. Then Il^f = Tlfcfi + l/c/2 is the orthogonal projection 



of H(n) to U{K), similarly to Q 
Our main result is the following. 



Theorem 2.1. 

particular, Id — 
of the form 

(2.10) 



Let K, satisfy (2.9). Then IL^Ai = Id- K in H D (K), with \\K\ 



H D (K) 



< 1. In 



K is invertible on Hrj(IC), and Ai restricted to Hrj(IC) has an explicit left inverse 



f=J2 ^njcAi/i, h = Ai/. 

m=0 



Remark 2.1. As discussed in the Introduction, K is not a compact operator as in the case of 
smooth sound speed. It follows from the proof of the theorem that the least upper bound of its 
essential spectrum (always less that 1) corresponds to the maximal portion of the high-frequency 
energy that is still held in Q at time t = T. 



Remark 2.2. Compared to [26, Theorem 1], there is a slight improvement in the theorem above; 
we reduce the lower bound of T by a factor of 2. Indeed, the proof also works in the case when c 
is smooth everywhere. In |26l Theorem 1], we assumed T > T(f2), where the latter is the length of 
the longest geodesic in Q. In particular, T(fl) is the diameter of f2, if the metric c~ 2 dx 2 is simple, 
see [25] • In the theorem above, T has to be larger than the "radius" of f2, see (2.8). Therefore, if 
T > T(Q)/2, for any (x,£) E 50, at least one of the geodesies 7^ and 7x,-£ would exit Cl at time 
T, hence (2.8) holds. Then K is still a contraction but not a compact operator anymore in general, 



even when c is smooth. Thus one can expect a slower convergence of (2.10) if T(0)/2 < T < T(Q). 



Remark 2.3. Consider the case now where K, does not satisfy (2.9). If there is an open set of 



singularities that does not reach <9f2, a stable recovery is impossible [24J. In either case however, 



a truncated version of the series (2.10) would provide an approximate parametrix that would 



recover the visible singularities, i.e., those in 14. By an approximate parametrix we mean a pseudo- 
differential operator elliptic in IA with a principal symbol converging to 1 in any compact in that 



set as the number of the terms in (2.10) increases. This shows that roughly speaking, if a recovery 



of the singularities is the primary goal, then only those in IA can be recovered in a "stable way" , 



and (2.10) works in that case as well, without the assumption (2.9). Details will be given in a 



forthcoming paper. 



3. Preliminaries 

Notice first that c 2 A is formally self-adjoint w.r.t. the measure c~ 2 dx. Given a domain U, and 
a function f = /2], define the energy 

E v (f)= [ {\V x fi\ 2 + c- 2 \f 2 \ 2 )dx. 
Ju 

Given a scalar function u(t,x), we also set u = [it,itt]; and then 

Eu(u(t))= [ (|V^| 2 + c- 2 |u t | 2 )dz. 
Ju 
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Here and below we use the notation u(t) = u(t, ■). We define the energy space %{U) as the 
completion of Cq°(U) x Cq°(U) under the energy norm 

||f f n(u) = E v {i) = J {\V x h\ 2 + c- 2 \h\ 2 ) dx. 

In particular, we define the space Hd(U) to be the completion of C^°(U) under the Dirichlet norm 

(3-1) \\f\\H D (u)= I |V^| 2 dx. 

J u 

The energy norm is topologically equivalent to the same norm with c = 1. It is easy to see 
that Hd(U) C H l {U), if U is bounded with smooth boundary, therefore, Hd{U) is topologically 
equivalent to Hq(U). If U = R n , this is true for n > 3 only [18 1 . By finite speed of propagation, 
the solution with compactly supported Cauchy data always stays in H 1 even when n = 2. This 
norm defines the energy space 

%{U)=H D {U)@L 2 {U). 
Here and below, L 2 (U) = L 2 (U; c~ 2 dx). Note also that 

(3-2) ||/||^ {[7) = (-c 2 A/,/) L2 . 

The wave equation then can be written down as the system 

(3-3) u, = Pu, P=( c2 ° A J), 

where u = [u, Ut] belongs to the energy space T~L. 

Proposition 3.1. Let U be an open subset o/R n with smooth boundary. Then the operator P with 
domain D(P) = {f E %{U); fa E Hd(U), Afa E L 2 (U)} is a unbounded skew- self adjoint operator 
on H(U) ifceL 00 , and c" 1 E L°°. 

Proof. Clearly, P is correctly defined on D(P), and the latter is dense in 7i(U). We will show first 
that P* d -P, i.e., that D(P*) D D(P) and P* = -P on £>(P). Let f, g be in £>(P). Then 

(3-4) (Pf , e)h{U) = ! (V/ 2 • V51 - V/i • Vg 2 ) dx. 

Ju 

This follows by integration by parts that can be easily justified by approximating g with smooth 
functions, see also [181 Theorem V.1.2]. Replace f by and f, and g by g to get 

(Pf,g)« (l0 = -(f,Pg)„ (t0 , Vf,g E D(P). 

We will show next that P* C —P. For that, in view of what we already proved, it is enough to 
show that D(P*) = D(P). By definition, g E D(P*) if and only if there exists h E H{U) so that 

(Pf,g)„ (t0 = (f,h) w(t0 , VfGD(P), 

see e.g., [231 §8-1]. This equality can also be written as 

/ (V/ 2 • V51 + (Afa)g 2 ) dx = f (VA • Vh + c- 2 fah 2 ) dx. 

Note that the l.h.s. is independent of c. Set h" = [hi,c~ 2 fi2\- It belongs to H. if and only if ( does. 
We therefore reduced the problem D(P*) = D(P) to its partial case when c = 1. In that case, this 
statement is well known, see e.g., |18} Theorem V.1.2]. In other words, such hJ E 7ii(U) exists, 
and equals —Pig, where "Hi, Pi correspond to c = 1, if and only if g E D(Pi) = D(P). Then we 
also get that h E H(U) exists (and, of course, equals — Pf). □ 
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Remark 3.1. If U is bounded then D(P) is topologically equivalent to H 2 (U) © If U is 

unbounded but suppf is a compact in U, then any f G D(P) belongs to that space as well. The 
case of non-compactly supported f is more delicate, see [18]. In this paper, we will deal with either 
U = R n or U = O; and in the former case, f is compactly supported. In the latter case, the 
definition of Hd(U) reflects Dirichlet boundary conditions on <9f2. 

By Stone's theorem, P is a generator of a strongly continuous group e* p of unitary operators on 
T~L{U) that leaves -D(P) invariant. For any f € D(P), we have g^e* p f = Pe tp f. In particular, for 
the first component u(t,x) of e* p f, we have that uu(t, •) is a continuous function of t with values 
in L 2 (U); the same holds for c 2 (-)Au(t, •), and uu = c 2 Au. 

Remark 3.2. An alternative proof of the proposition is to use well-posedness results in [32], for 
example, to show that the solution is given by a strongly continuous unitary group; and deduce 
from there that its generator must be skew-selfadjoint. 

For f G D(P), the following transmission conditions hold for f\ 



( 3 - 5 ) ^lr int =/i|r cxt < 



dfx df 
d 



v 



Vint 3v 



that just reflects the fact that f\ is in H locally. 



4. Geometric Optics 

4.1. A parametrix of the Cauchy problem with data at t = 0. We start with a standard 
geometric optics construction. We assume first that c is smooth in the region where we construct 
the parametrix. Fix (xq,^°) £ T*Q \ 0. We identify vectors and covectors away from T by the 
metric c~ 2 dx 2 . Let u be a solution of the wave equation with initial data [it, ut] = [/i,/2], with 
f = [fi, fj[ having a wave front set in some small conic neighborhood of (xo,£ )- 



In a neighborhood of (0,xq), the solution to (2.1) is given by 



(4.1) u(t,x) = (2vr)-" [ elMt ' X '° (*i*M,t)M$ + a2Ax,Z,t)M0) d ^ 

modulo smooth terms. Here the phase functions <f>± are positively homogeneous of order 1 in £ and 
solve the eikonal equations 

(4.2) ±d t <f>± + c(x)\V x (f>±\ =0, cp±\ t=0 = x-Z, 

while a,j ± are classical amplitudes of order solving the corresponding transport equations along 
bicharacteristics issued from some conic neighborhood of (x,±£) for all (x,£) € WF(f), see (H 
p. 128] or [311 eqn. (VI. 1.50)]. In particular, aj t ± satisfy 

ai )+ + a\- = 1, 02,+ + a.2- =0 for t = 0. 

Since dt4>± = = F C ( 2; )|£| f° r t = 0, and ut = fa for t = 0, we also see that for the principal terms 

of aj t ± ~ Y,k>o a j!± fc ^ we ^ ave 

4% = 4% c(x)iei(4°l-4°|) = i fart = o. 

Therefore, 

(4.3) afi = 4!l = lA 4°| = -4°1 = ^ M fort = 0. 
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Note that if c = 1, then (fr± = x ■ £ =F £|£|, and ai >+ = a\- = 1/2, a2,+ = —a-z - = (2|£|) 1 . The 
principal terms af\. satisfy the homogeneous transport equations 

(4.4) {{dt<t>±)d t - c 2 V^± • V x + C±) = 

with initial conditions given by ( |4.3[ ), where 

2C ± = (d 2 -c 2 A)«/> ± , 

see also gU eqn. (VI.1.49)]. 

By the stationary phase method, for a = +, singularities starting from (x,^) £ WF(f) propagate 
along geodesies in the phase space issued from (x,£), i.e., they stay on the curve (7x,§(i)> 7z,?W); 
and from (x, — £), for a = — , i.e., they stay on the curve (7 Xi _^(£), 7 Xi _^(t)). This is consistent with 
the general propagation of singularities theory for the wave equation because the principal symbol 
— t 2 + c 2 |£| 2 of the wave operator has two roots r = ±c|£|. 

The construction is valid as long as the eikonal equations are solvable, i.e., for \t\ small enough, 



at least. We call the corresponding parametrices in (4.1) u + and U-. It is well known that u + + u 



represents the solution of (2.1) up to a smoothing operator applied to f. 



4.2. Projections to the positive and the negative wave speeds. The zeros of the principal 
symbol of the wave operator, in regions where c is smooth, are given by r = ±c(:r)|£|, that we 
call wave speeds. We constructed above parametrices u± for the corresponding solutions. We will 
present here a functional analysis point of view that allows us to project the initial data f to data 
II-|-f, so that, up to smoothing operators, u± corresponds to initial data II-i-f. 

Assume that c(x) is extended form the maximal connected component of R n \ T containing xq 
to the whole R n in a smooth way so that < 1/C < c(x) < C. Let 

(4.5) Q = (-c 2 A) 1//2 , 

where the operator in the parentheses is the natural self-adjoint extension of — c 2 A to L 2 (R™, c~ 2 dx), 
and the square root exists by the functional calculus. Moreover, Q is an elliptic ^DO of order 1 in 
any open set; and let Q _1 denote a fixed parametrix. 



It is well known that the solution to (2.1) can be written as 



(4.6) u = cos(tQ)/! + Si ^/ 2 , 

and the latter operator is defined by the functional calculus as 4>(t, Q) with cp(t, A) = sin(i-)/- € C°°. 
Based on that, we can write 

(4.7) e tp = e itQ n+ + e^II., 
where 

It is straightforward to see that II± are orthogonal projections in %, up to errors of smoothing 
type. Then given f G % supported on Q, one has u± = e' p f±, with i± := II-i-f. 
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4.3. Analysis at the boundary. We will analyze what happens when the geodesic (xo, £°) issued 
from (xo,£°), xq g" T, hits T for first time, under some assumptions. Let the open sets Oi nt , $7 ex t ; 
be the "interior" and the "exterior" part of £1 near xq, according to the orientation of T. They only 
need to be defined near xq. Let us assume that this geodesic hits V from f2; nt . We will construct 
here a microlocal representation of the reflected and the transmitted waves near the boundary. 

Extend c\n int in a smooth way in a small neighborhood on the other side of T, and let u+ be 
the solution described above, defined in some neighborhood of that geodesic segment. Since we 
are only going to use u+ in the microlocal construction described below, and we will need only 
the trace of u + on R + x Y near the first contact of the bicharacteristic from (xo,£°) with T, the 
particular extension of c would not affect the microlocal expansion but may affect the smoothing 
part. 

Set 

(4.9) h:=u + \ Kxr . 

Let (ti,xi) G R+ x T be the point where the geodesic from j xo ^o hits T for the first time, see 
Figure [2} We assume that such t\ exists. Let £ x be the tangent covector to that geodesic at 
(ti,xi). Assume that ^° is unit covector in the metric c _2 dx, then so is £ x (in the metric c^ 2 dx), 
i.e., Ci n t|£| = 1, where |£| is the Euclidean norm. Assume that is transversal to V. In view of 



condition (2.9), this is the case that we need to study. 

Standard microlocal arguments show, see |26| Proposition 3] for details, that the map [fx, f<z\ 1— > h 
is an elliptic Fourier Integral Operator (FIO) with a canonical relation that is locally a canonical 
graph described in [261 Proposition 3]. That diffeomorphism maps (xo,£°) into (tx,xx, 1, (£ )'), 
where the prime stands for the tangential projection onto T*T; and that maps extends as a positively 
homogeneous one of order one w.r.t. the dual variable. In particular, the dual variable t to t stays 
positive. In fact, WF(«) is in the characteristic set t 2 — c 2 (x)|£| 2 = 0, and (x, £) belongs to some 
small neighborhood of (x±, t; 1 ). The wave front set WF(/i) is given by (x, £') G T*T, (x, £) G WF(u), 
where £' is the tangential projection of £ to the boundary. Then (t, x, r, £') is the image of some 
(x, £) close to (xo, £°) under the canonical map above. Here (x, ^) is such that the x-projection x(s) 
of the bicharacteristic from it hits T for the first time at time for the value of s given by sc{x) = t. 
Since r 2 - c 2 nt (x)|£| 2 = 0, for the projection £' we have r 2 - c 2 nt (x)|^'| 2 > 0, where (x,£') G T*T, 
and is the norm of the covector in the metric on T induced by the Euclidean one. 

The microlocal regions of T*(R xT) 3 (t,x, r, £') with respect to the sound speed c; nt , i.e., in 
Cli Q t, are defined as follows: 

hyperbolic region: Ci nt (x)|^'| < r, 
glancing manifold: Ci n t(x)|£'| = r, 
elliptic region: Ci nt (x)|£'| > r. 

One has a similar classification of T*T with respect to the sound speed c ex t- A ray that hits T 
transversely, coming from Qi n t, has a tangential projection on T*(R x T) in the hyperbolic region 
relative to c; n t- If Qnt < c ox t, that projection may belong to any of the three microlocal regions 
w.r.t. the speed Cj n t. If c- m t > c ex t, then that projection is always in the hyperbolic region for c ex t • 
When we have a ray that hits T from f2 e xt> then those two cases are reversed. 

4.4. The reflected and the transmitted waves. We will analyze the case where belongs to 
the hyperbolic region with respect to both q n t and c ex t, i-e., we will work with £' in a neighborhood 
of (t; 1 )' satisfying 

(4-10) c^r 2 - |e'| 2 > 0, c- 2 r 2 - |e'| 2 > 0. 
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The analysis also applies to the case of a ray coming from fi e xtj under the same assumption. We will 
confirm below in this setting the well known fact that under that condition, such a ray splits into a 
reflected ray with the same tangential component of the velocity that returns to the interior S7i nt , 
and a transmitted one, again with the same tangential component of the velocity, that propagates in 
^ext ■ We will also compute the amplitudes and the energy at high frequencies of the corresponding 
asymptotic solutions. 

Choose local coordinates on T that we denote by x 1 , and a normal coordinate x n to T so that 
x n > in f^ext, and \x n \ is the Euclidean distance to T; then x = (x',x n ). We will express the 
solution u + in R x Oi n t that we defined above, as well as a reflected solution ur in the same set; 
and a transmitted one ut in R x ex t, up to smoothing terms in the form 

u a = (27r)-" f jO'^rt'Kit, x, r, £')h(T, O dr de', a = +, R, T, 



(4.11) u c 

where h := ,f RxR n-i e~ 1 ^~ tTJ> " x ''^'^h{t, x')dt dx' . We chose to alter the sign of r so that if c = 1, then 
the phase function in (4.11 ) would equal (p+, i.e., then ip + = —tr + x • £. The three phase functions 



cp+, ipn, ifT solve the eikonal equation 
(4.12) dt<p a + c(x)\X7 x <p a \ 



0, <fo 



-tr + x' ■ f. 



The right choice of the sign in front of dtip+, see (4.2), is the positive one because dt<p+ = —t < 
for x n = 0, and that derivative must remain negative near the boundary as well. We see below that 
<PR t T have the same boundary values on x n = 0, therefore they satisfy the same eikonal equation, 
with the same choice of the sign. 

Let now h be a compactly supported distribution onRxT with WF(/i) in a small conic neighbor- 
hood of (t\,x\, 1, (C 1 )'). We will take h as in (4.9) eventually, with u + the solution corresponding 
to initial data f at t = but in what follows, h is arbitrary as long as WF(/i) has that property, 
and u+ is determined through h. We now look for a parametrix 



(4.13) 



U + + UR + UT 



near (t\,xi) with u+, ur, ut of the type (4.11), satisfying the wave equation and (4.9). We use the 
notation for u + now for a parametrix in f2i nt having singularities that come from the past and hit 
r; i.e., for an outgoing solution. The subscript + is there to remind us that this is related to the 
positive sound speed c(x)|£|. Next, ur is a solution with singularities that are obtained form those 
of u+ by reflection; they propagate back to fimt- It is an outgoing solution in Sl; nt . And finally, ut 
is a solution in Q ext with singularities that go away from T as time increases; hence it is outgoing 



there. To satisfy the first transmission condition in (3.5), we need to have 
(4.14) ip T = ip R = ip + = -tr + x ■ £' for x n = 0, 



that explains the same boundary condition in (4.12), and 



(4.15) l + b R = b T forx n = 0. 

In particular, for the leading terms of the amplitudes we get 



(4.16) 



1 for x n = 0. 



To satisfy the second transmission condition, we require 



(4.17) 
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Expanding this in a series of homogeneous in (r, £) terms, we get series of initial conditions for the 
transport equations that follow. Comparing the leading order terms only, we get 



(4.18) 



d <PT b (0) 

dx n T 



d( PR b (Q) 

Q x n R 



for x r ' 



0. 



The linear system (4.16), (4.18) for b R 
(4.19) 



(0), 



*=0, 

difT 
dx n 



dx r 

=o has determinant 

9<pr\ 
1 I 

/ x n =0 

i,( )l 



0. 



Provided that this determinant is non-zero near x\, we can solve for b R '\ x n=o, bj! \ x n=o- Moreover, 

f— ?)i (—1)1 

the determination of each subsequent term b R '\ x n =0: bj, \x n =o m the asymptotic expansion of 
x n =o, br\x n =o can be found by (4.17) by solving a linear system with the same (non-zero) 
determinant. 



4.5. Solving the eikonal equations. As it is well known, the eikonal equation (4.12) in any 
fixed side of R x T, near (ti,xi), has two solutions. They are determined by a choice of the sign 
of the normal derivative on R x T and the boundary condition. We will make the choice of the 
signs according to the desired properties for the singularities of u+, ur, ut- Let V x ' denote the 
tangential gradient on T. By (4.14), 

(4.20) 



V., 



■<PT 



Vx'VR = Vx><P+ = £' 3 dt<PT = d t <fR = d t cp + = -t for x n = 0. 



Using the eikonal equation (4.12) and the boundary condition there, we get 

dip + 



(4.21) 



C int r2 



le 



l\2 



for x r ' 



0. 



dt dx n 

We made a sign choice for the square root here based on the required property of u+ described 
above. This shows in particular, that the map h \— > du+/dt (that is just d/dt), and the interior 
incoming Dirichlet to Neumann map 

du + 



int,in 



: h H> 



dv 



Rxr 



are locally ^DOs of order 1 with principal symbols given by 



At, and 



(4.22) 



<7p(-Nint,in) 



dx r 



1^ 



/|2 



The notion "interior incoming" is related to the fact that locally, near (ii,xi), we are solving a 
mixed problem in R x fij nt with lateral boundary value h and zero Cauchy data for t > 0. 

Consider ip r next. The reflected phase <£r solves the same eikonal equation, with the same 
boundary condition, as <p+. By the eikonal equation (4.12), we must have 

d<PR 



(4.23) 



dx T - 



9<P-\ 
' dx Ti 



for x n = 0. 



The "+" choice will give us the solution (p + for ipR. We chose the negative sign, that uniquely 



determines a solution locally, that we call ipR, i.e., 
(4.24) 9iPR - ^ + 



for x n = 0. 



dx n dx r ' 

Therefore, Vy^ on the boundary is obtained from V</3+ by inverting the sign of the normal deriv- 
ative. This corresponds to the usual law of reflection. Therefore, 

dif R _ difR 



(4.25) 



dt 



dx r - 



le 



/|2 



for x r 



0. 
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Figure 2. 

In particular, 8ur/ 'cte ra |R X r can be obtained from ii,R|Rxr> that we still need to determine, via the 
interior outgoing Dirichlet-to-Neumann map 

du R 

Rxr dx n Rxr 

that is locally a first order 'I/DO with principal symbol 



^Vint,out : U R 



(4.26) 



0"p(^V"int,out) = 1 



■ d<PR 



dt 



V^ 2 - Id 2 - 



To construct ifT, we work in 0, ex t- We define ipx as the solution of (4.12) with the following 



choice of a normal derivative. This time ipx and 93+ solve the eikonal equation at different sides of 
r, and c has a jump at F. By ( 4.20| ), 



(4.27) 



-'ext 



IC'| 2 + 



dtpi 



dx r 



t 2 for x n = 0. 



We solve this equation for \difT/dx n \ 2 . Under the assumption (4.10), this solution is positive, 
therefore we can solve for dipx/dx™ to get 

dipT 



(4.28) 



dx n 



c- 2 t r2 _| f |2 forx - = . 



The positive sign of the square root is determined by the requirement the singularity to be outgoing. 
In particular, we get that the exterior outgoing Dirichlet-to Neumann map 

8ut 

Rxr 



dx r 



Rxr 



-Next, out : ut 

has principal symbol 
(4-29) a p (iV ext , out ) = i|^ = Vc c - 2 t r2_| e 



r\2 
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For future reference, we note that the following inequality holds 

Qnt 



(4.30) 



d<p T d<p 
< tt-t < 7 



dx r 



dx 11 ' 



7 : = 



max < 1 . 

r c ext 



4.6. Amplitude and Energy Calculations. By (4.25), (4.28), the determinant (4.19) is nega- 
tive. Solving (4.16) and (4.18) then yields 

(4.31) 



1,(0) 



2dp + /dx n 



dip + 1 dx n + difT I dx n ' 



dip + /dx n — dipr/dx™ 
dip + 1 dx n + difT I dx n 



for x n = 0. 



As explained below (4.19), we can get initial conditions for the subsequent transport equations, 
and then solve those transport equation. By (4.14), the maps 

(4.32) Pr- u R \n.xr, Pt ■ h mtIrxT 



are ^DOs of order with principal symbols equal to , 6j? restricted to R x T, see (4.31). We 
recall (4.9) as well. 

We estimate next the amount of energy that is transmitted in f2 e xt- We will do it only based 
on the principal term in our parametrix. That corresponds to an estimate of the solution operator 
corresponding to tra nsmi ssion, up to com pact operators, as we show below. 

A quick look at (4.31), see also (4.16) shows that bj! > 1- This may look strange because 
we should have only a fraction of the energy transmitted, and the rest is reflected. There is no 
contradiction however because the energy is not proportional to the amplitude. 

Let u solve [d^ — c 2 A)n = in the bounded domain U with smooth boundary for tf < t < t" 
with some t' < t". A direct calculation yields 

(4.33) Ev{u{t")) = Eu(u(t')) + 2M [ <>U 

J[t',t"]xdU 

We will use this to estimate the energy of ut in f2 ex t- Since the wave front set of ut is contained 
in some small neighborhood of the transmitted bicharacteristic, we have smooth data for t = 0. 
Therefore, if t2 > h is fixed closed enough to t\, we can apply (4.33) to a large ball minus Sl; n t to 
get that modulo a compact operator applied to h, 



ut— dt dS. 
ov 



(4.34) 

Therefore, 
(4.35) 



^ cxt (u T (t 2 )) = 2K 



[o,t 2 ]xr 



But Out 
dt dv 



dtdS. 



En ext (u T (t 2 )) = 2R(P t u T ,N extput u T ) = R(2P^N: xtout P t P T h,h), 



where (•, •) is the inner product in R x R n , and Pt = d/dt. 

Apply similar arguments to u + in Oj nt . Since the bicharacteristics leave f2; n t) we have modulo 
smoother terms 



(4.36) 



0-^ int (u + (0)) + 25R 



[o,t 2 ]xr 



dt dv 



dtdS. 



Similarly we get, see again (4.32), 

(4.37) £n int (u + (0)) = -2R(P t h, N- mt - m h) = U(2N^ t>hl P t h, h). 

For the principal symbols of the operators in (4.35), (4.37) we have 
a p (2P*N*. 



(4.38) 



tPtPr) dipx/dv 



° P m nt . m p t ) 



dip + /dv 



,(°) 



4(d(p + / dh>)(d(fT / dv) 
(d(p + /dv + difx/dv) 2 
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Denote for a moment a := dip^/dv, b := d(px/dv. Then the quotient above equals Aab/(a + b) 2 < 1 



that confirms that the reflected wave has less energy than the incident one. By (4.30), < b < 7a, 
< 7 < 1. This easily implies 

4a6 47 

Therefore, the expression in the middle represents an upper bound of the portion of the total energy 
that gets transmitted in the asymptotic regime when the frequency tends to infinity. To get a lower 
bound, assume in addition that b > bo > and a < ao for some ao, bo, i.e., 

(4.40) 0<6 <^, ^±<a . 



dv dv 



Then 



4ab 46^/7 

4.41 - > - 0/ i > 0. 

(a + 6) 2 (l + 7 ) 2 a 2 , 

This is a lower bound of the ratio of the high frequency energy that is transmitted. As we can see, 
if the transmitted ray gets very close to a tangent one, that ratio tends to 0. 

So far this is still not a proof of such a statement but just a heuristic argument. We will formulate 



and prove this below in Proposition 5.1 



4.7. Snell's Law. Assume now that is in the hyperbolic region for q n t but not necessarily 
for Cext- This corresponds to a ray hitting V from the "interior" Oj nt . There is no change in solving 
the eikonal equation for ipn but a real phase <px does not exist if the expression under the square 



root in (4.28) is negative. This happens when is in the elliptic region for c ex t- Then there is 
no transmitted singularity in the parametrix. We analyze this case below. If c; n t > c ex t, then (£ )' 
that is in the hyperbolic region for c; n t by assumption, also falls into the hyperbolic region for the 
speed c ext , i-e., there is always a transmitted ray. If ci nt < c ex t, then existence of a transmitted 
wave depends on where belongs w.r.t. c ex t- 

Let a be the angle that £ = dip+/dx n makes with the (co)-normal represented by dx n , and let 
f3 be the angle between the latter and £t '■= dipx/dx™ . We have 

(4.42) = le 1 ! sina = ctJt sin a, |£'| = |£ T | sin/3 = c^rsin/3 



By (4.42), we recover Snell's law 

sin a Ci n t 
sin/3 c ext ' 



(4.43) 



Assume now that Ci n t < c ex t, which is the case where there might be no transmitted ray. Denote 
by 

(4.44) ao(x) = arcsin(ci n t(x)/c cx t(x)) 

the critical angle at any x S F that places (i; 1 )' in the glancing manifold w.r.t. c ex t- Then the 
transmitted wave does not exist when a > ao', more precisely we do not have a real phase function 
ipT in that case. It exists, when a < ao- In the critical case a = ao, this construction provides an 
outgoing ray tangent to T that we are not going to analyze. 
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4.8. The full internal reflection case. Assume now that (t; 1 )' is in the elliptic region w.r.t. 
Cextj then there is no transmitted singularity, but one can still construct a parametrix for the 
"evanescent" wave in O ext ; and there is a reflected ray. This is known as a full internal reflection. 
We give details below. 

We proceed as above with one essential difference. There is no real valued solution fT to the 



eikonal equation (4.12) outside Qq- Similarly to (4.28), we get formally. 



(4.45) 9 -^= ^\e\ 2 -c-^ forx" = 0. 

The choice of the sign of the square root is dictated by the requirement that the so parametrix 
( |4.11 ) with g = T be exponentially decreasing away from F instead of exponentially increasing. 



In general, the eikonal equation may not be solvable but one can still construct solutions modulo 
0((x n )°°). The same applies to the transport equations. One can show that the 0((x n )°°) error 



does not change the properties of ut to be a parametrix. In particular, in (4.35) in this case one 
gets 

(4-46) E Qcxt (u T (t 2 )) = 0, 



because the principal term of dux/dv in (4.34) now is pure imaginary instead of being real. More- 
over, ut is smooth in f2 e xt- Therefore, no energy, as far as the principal part only is considered, is 
transmitted to J7 e xt- That does not mean that the solution vanishes there, of course. 

4.9. Glancing, gliding rays and other cases. We do not analyze the cases where is in the 
glancing manifold w.r.t. to one of the speeds. We can do that because the analysis of those cases in 
not needed because of our assumptions guaranteeing no tangent rays. The analysis there is more 
delicate, and we refer to [291 123 EI] for more details and examples. We do not analyze either the 
case where {^})' is in the elliptic region with respect to either speed. 

4.10. Justification of the parametrix. Denote by ur = [uR,dtUn\, uy = [uT,dtUT] the approx- 
imate solutions constructed above, defined for t in some neighborhood of t 2 . Then u^j = Vr/j, 



ut = Vy/i, where Vrt are the FIOs constructed above. Let u + be the solution of (2.1) defined 
above, with initial data II + f at t = having wave front set in a small neighborhood of (xo,^°). 
The map A + : f i— > u+IrxT = h is an FIO described in [26]. Then near (ti, xi), 

u R = V R Af , u r = VjAf , 

the former supported in R x Oi n t, and the later in R x f2 C xt- So far we had two objects that we 



denoted by u+: first, the parametrix of the solution of (2.1) corresponding to the positive sound 
speed c(x)|£|; and the parametrix in R x f2j nt for the incoming solution corresponding to boundary 
value h. When h = A + f, those two parametrices coincide up to a smooth term, as it is not hard 
to see (the second one is a back- projection and is discussed in |26j . in fact). This justifies the same 
notation for them that we will keep. 

Consider the parametrix v p := u+ + ur + ut- We can always assume that its support is in some 
small neighborhood of the geodesic that hits R x T at (ti,xi) and is tangent to ^ 1 there; and then 
reflects, and another branch refracts, see Figure [2] In particular, then h has f-support near t = ti, 
let us say that this included in the interval [ti — e, ti + e] with some e > 0. At t = ti := t\ + 2e, let 
X2 be the position of the reflected ray, and let £ 2 be its unit co-direction. Then WF(n^j(t2, •)) * s 111 
a small conic neighborhood of (x2,£ 2 ). 

Let v(t, ■) = e* p II_|_f be the exact solution, with some fixed choice of the parametrix Q~ 1 in the 
definition of n + , properly supported. Consider w = v — v p in [0, £2] x R n - It satisfies 
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(4.47) (df - c 2 A)w\ [0MxClint e (d 2 t - c 2 AH [CM2]x ^ xt G 

dw I dw I 

(4.48) ^l[o,t 2 ]xr cxt -^l[o,t 2 ]xr int 6 C°°, ^ l[o,t 2 ]xr cxt " ^l[o )i2 ]xr int e C °°- 

On the other hand, for < t <C 1, v is smooth. Let x G C°°(R) be a function that vanishes in 



(—oo,5] and equals 1 on [25, oo), < (J < 1. Then := x(t)w(t,x) still satisfies (4.47), (4.48) and 
also vanishes for t < 0. By [32j Theorem 1.36], w is smooth in [0,t2] X ^int> up to the boundary, 
and is also smooth in [0,i2] x ^ext; up to the boundary. Therefore, 

(4.49) v(i,-) = v P (t J -)+K t f, 

for any t G [0, £2]) where Kf is a compact operator in T~L, depending smoothly on t. The operator Kf 
depends on Q as well. Therefore, the parametrix coincides with the exact solution up to a compact 
operator that is also smoothing in the sense described above. 

5. Proof of the main result 

We start with a lemma, that in principle is well known and is related to the classical Dirichlet 
principle. 

Lemma 5.1. Let Q be a bounded set with smooth boundary. Define the map 

as follows: Hq u — xl — (pu? where (p u is the solution to 

Acf> u = in ft, (f) u \an = u\dn- 

Then 

H n fi"llH D (fi) < \\ u \\H D (n), 
where \\u\\h d (q) is the Dirichlet norm of u extended to functions in H 1 ^) that may not vanish on 

Set II[)U = [Uqui,u 2 ]. Then 

E n (n n u) < E n (u). 



Proof. In what follows, (•, -)H D (n) is the inner product in Hr>(£l), see (3.1), applied to functions 
that belong to i? 1 (Q) but maybe not to Hr>(Q) because they may not vanish on dQ,. Since u = (\> u 
on dil, and since 4> u is harmonic, we get 



(u - (j) u , (pu)H D (n) = I V(n - (j) u ) ■ V(f) u dx = 0. 

Then 



n 



\\ u - 4>u\\H D (n) — \\ u \\H D (n) - W&uWhoW ^ 

The second statement of the lemma follows immediately from the first one. □ 

It is easy to see that Hq is a projection, orthogonal with respect to the product (•, -)H D (n)- This 
products defines a seminorm on i7 1 (0) only. 



Before giving the proof of the main theorem, we will show that condition (2.9) implies the 
following local energy decay estimate. 
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Proposition 5.1. Let (2.9) be satisfied. Then there exists < [i < 1 so that 
(5.1) E n (u(T)) < fiE n ({) 

for any solution u with initial data f = [/i,0] E 'H(IC). 

Proof. We will use the geometric optics construction in section |4j Let X = diag(X, X) be a zeroth 
order ^DO with small enough essential support supported near some (xo, £°) E T*K. \ 0. Let 
v := e* p Il + Xf be the corresponding solution. Recall that n + restricts propagation of singularities 
to the positive wave speed c(x)|£| only, while X localizes near (xo,£ )• Assume for now, that the 
ray through (xq,£ ) gives rise to both a reflected and a transmitted one, as in section 



4.4 



We will use the energy computation in (4.35) and (4.37), under the assumption (4.40). According 
to those relations, (4.38) and (4.41), 

(5-2) E n ^ t (u T (t 2 )) - ^ int (u + (0)) = »(M/i, h), 

with M a ^DO of order with non-negative principal symbol, some ^ E (0,1), and we use the 
notation in that section. By the Garding inequality, 

(5.3) £n cxt (uT(i 2 )) > //#n int (u + (0)) - C\\hf Hl/2{Rxdn) = ^'E n . int (H + Xf ) - C\\hf Hl/2{Rxdny 

The map A.i + : [fi, 0] i— > h is an FIO of order with a canonical relation of graph type discussed 
above. Similarly, the map A2 5 + : [0, f 2 ] h > h is an FIO of order —1 with the same canonical relation. 
Then 



\h\ 



H 1 / 2 {R,xdVt) 



< C\\u, 



Therefore, 

(5-4) ^Wt(ur(*2)) > ^'^(n+xf) - cun+xf f m/ ^ H - 1/2 . 

Next, assume that the ray through (xo,^°) gives rise to a reflected ray only, as in section 
Then by fl4.46| ), 

(5.5) 



4.8 



£;n ext (u T (i 2 ))<c||n + xf||^ 1/2eH . 



1/2- 



By energy preservation, 
(5-6) E nc 



: (u fl (t 2 ))>£?n iIlt (n + Xf)-C||n. 



-Xf ||^l/2 ( 



and we used the fact here that u + (t2,-) can be obtain from the initial conditions by applying a 
smoothing operator. 

We will apply those arguments now in a more general situation. Assume now that (x ,£°) E U, 
see (2.8). By the definition of U, there is a path starting from (xo,£°) or from (xo,^°) so that 
it consists of finitely many geodesic segments; and at least one of them is a geodesic hitting T 
transversely for the first time at some x\ an angle strictly greater than the critical one tt/2 — ao at 
x\. Without loss of generality, we can assume that we have the former case: the path starts from 
( x o>£ )- Then the next geodesic segment in that path is a transmitted one that is not tangent to 
r. The inequality (4.40) then holds for the corresponding phase with some bo > 0. 



Assume that the essential support of the symbol q is so small that the analysis in section 4.8 



applies near each full internal reflection. Then step by step, applying (5.6) consecutively, we can 
reduce the analysis to the case where x\ is the first point where the geodesic from (xq,^°) hits T. 
Then we apply (5.4) to estimate the amount of energy that has been transmitted on the other side 
of r. Let for a moment assume that the transmitted ray leaves fi for time t < T without further 
contact with T. Then for any t > t±, not necessarily close to it, the map ux(ii,-) \-t ux(t, •) is 
an FIO that represents the solution of the wave equation outside Oo corresponding to the positive 
wave speed. In particular, it is not affected by the presence of a jump at T. Therefore, an estimate 
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equivalent to (5.4) is preserved for such t as well but now the energy of u(T) is concentrated outside 
fi, up to smoothing terms. Let fio be an open set with a smooth boundary so that 1C C fio <s fi\r. 
Then 



(5.7) 



ixn+fl 



\n(u ) - CE ~R- n \n( u { T )) + C||f 

Above, we also modify II + and X by smoothing operators, if needed, so that XII + f belongs to the 
indicated energy space, i.e., they vanish in R n \ Qq. 

Consider now the case where the ray guaranteed by (2.8), (2.9) may hit T again, even more than 
once, before leaving fi and not coming back. At each such event, there will be no loss of energy 
at high frequencies, as in (|5.5[ ), or there will be positive portion of the transmitted high frequency 
energy, as in (5.4). Then (5.7) is still true. 

So far we did not use the assumption /2 = 0. Now, since f = 0], we get n + Xf = iQXfx], 
see (4.8). By the ellipticity of Q, 

( 5 - 8 ) \\ x fi\\H D (Q ) < C f ||u(^))||fl-i(R»\n)©i2(a«\n)) + C|I/iIIhi/2- 

By a compactness argument, in a conical neighborhood of S*tC, we can take a finite pseudo- 
differential partition of unity 1 = Y2j Xj °f symbols of ^DOs Xj localizing in conical neighborhoods 
of a finite number of points (xj,^) G S*JC. Thus we get 

( 5 - 9 ) \\h\\H D (n ) < C|l u ( r ))ll// 1 (R n \n)eL 2 (R n \n)) + c \\h Wh 1 / 2 ■ 

Consider the bounded map 

(5.10) H D {K) 3 fi i— »• u(T) G H\H n \ fi) © L 2 (R n \ fi)). 

We claim that it is injective. Indeed, assume that for some f G /C, for the corresponding u we have 

u{T, x) = 0, for x £ Q. 
By finite domain of dependence in R x R" \ Q, where c = 1, we get 

(5.11) u(t,x) = when dist e (x, O) > \T — t\, 

where dist e stands for the Euclidean distance. One the other hand, we also have 

(5.12) u(t, x) = when dist e (x,f2) > \t\. 

Here, we applied finite domain of dependence argument outside fi, as well. Note that this is not 
sharp, because at least when c does not jump at T, then the Euclidean distance dist e (x, Q) can be 
replaced by the distance in the metric c~ 2 dx 2 between /C and x. 
Therefore, 

(5.13) u(t, x) = when dist e (x, Q) > T/2, -T/2 < t < 3T/2. 

The solution u extends to an even function of t that is still a solution of the wave equation because 
/2 = 0. Then one gets that (|5.13 ) actually holds for \t\ < 3T/2. We will conclude next by the 
unique continuation Theorem 6.1 that u = on [0, T] x fi, therefore, fx = 0. 

To this end, notice fist that from John's theorem (equivalent to Tataru's unique continuation 
result |26} Theorem 2] in the Euclidean setting), we get u = on [—T,T] x R n \ fi. Fix xq G fi. 
Then there is a piecewise smooth curve starting at xq in direction either £° or — £°, where £° is 
arbitrary and fixed, of length less than T that reaches 3fi. This means that dist(xo, R n \ ^) < T, 
see Definition 6.1 Then by Theorem 6.1 u(0, •) = near xq. Since xq was arbitrary, we get / = 0. 
This completes the proof of the claim that (5.10) is injective. 
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t=3T/2 



t=T/2 



Figure 3. 

Now, by [321 Proposition 5.3.1], since the inclusion He>(IC) H 1 / 2 is compact, there is an 
estimate as in (5.9), with a different C, with the last term missing, i.e., 

( 5 - 14 ) \\fl\\H D (K.) < C|l u ( T )llHi(R™\^)eL 2 (R' i \n)- 

By finite speed of propagation, u(T) is supported in some large ball B. Apply the Poincare 
inequality [5] to conclude that 

(5-15) \\h\\H D{ K) < CE Rn \ n (u(T)). 



Then, using Lemma 5.1 we get 

E n (u(T)) = £7r»(u(T)) - E Rn \ n (u(T)) 
< E K (f) - E K (f)/C 
= (1 - 1/C 2 )^(f), 

and this completes the proof. 



□ 



Proof of Theorem 2A_ Let f G , H{Q) first. Let w solve 

f (5 4 2 -c 2 A)u7 = 



(5.16) 



in (0,T) x n, 



w \[0,T]xdQ. 
w\t=T 
W t \t=T 



0, 

«|t=T " 



where u solves (2.1) with a given f EH, and (/> is as in (2.4) with /i = A/. We can also write 
(5-17) [w,w t }\t=T = n c ([u,ut\\{ T }xn) , 



where IIq : H 1 ^) — )■ is the projection introduced in Lemma 5.1 Then 

(5.18) E n (w(T)) < E a (u(T)). 



Let u be the solution of (2.4) with /i = A/. Then t> + u; solves the same initial boundary value 



problem in [0,T] x 17 that u does (with initial conditions at t = T), therefore u = v + w. Let now 
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f = [fx, 0] G %{K). Restrict u to t = and project to Hd(JC) to get 

/i = n^iAi/i + n^(o,o. 

Set 

K/i = njc«;(0, •)• 

It remains to show now that K extends to a contraction. 



By (5.18) and Proposition 5.1, the energy of the initial conditions in (5.16) satisfies the inequality 
(5-19) £ n (w(T))< Mil/if^. 

Since the Dirichlet boundary condition is energy preserving, we get 
(5-20) E u (w(0)) = E n (w(T)) < Hi/ill^- 

Therefore, 

(5.21) \\Kfx\\ 2 HD{)C) < E K (n K w(0)) < E K (w(0)) < £ n (w(0)) < n\\fi\\ 2 Hj3 . 

This completes the proof. □ 

6. Unique continuation 

Definition 6.1. Given x, y in R™, let dist(x,y) be the infimum of the length of all smooth curves 
connecting x and y, intersecting V transversally at each common point, where the length is measured 
in the metric c~ 2 dx 2 . 

The next result is a unique continuation theorem. The proof is based on the smooth case. 
Theorem 6.1. Let xq g" £1 and T > 0. Assume that [u,ut] G C(R; H\ oc ) and u satisfies 

(d 2 - c 2 A)u = 0, 

in a neighborhood of the set 

(6.1) U = {(t, x); \t\ + dist(x , x) < T} 
and vanishes in a neighborhood of[—T,T] x {xq}. Then 

(6.2) u(t,x) = inU. 

Proof. When c is smooth, the theorem follows from Tataru's unique continuation result [27, 28j, 
see |26] . In the case we consider, we will base our proof on the smooth case. 

We show next that one can assume without loss of generality that [u,ut] G C(R; D(P)), where 
D(P) stands for the domain of P equipped with the graph topology. Assume first that u solves the 
wave equation in the whole R x R™. By finite speed of propagation [32], we can always assume 
that u has compactly supported Cauchy data for t = 0. Then u = e tp f, with some f € %■ Take 
4> £ C^°(R) and consider the convolution = u *t <p °f u with <j) in the t variable. Then is a 
solution of the wave equation with support close enough to that of u, if supp <fi is small enough; so 
if we prove the theorem for u^, that would be enough. The initial value for u<^ = [u^, dtu^] at t = 
is given by u^(0, •) = f u^(— s, -)(f>(s) ds. Apply the generator P to that to get 

Pu (O,-) = J Pu (-s,>(s)ds = J u 4> (-s,-)(f)'(s)dseH. 

Therefore, the initial condition now is in the domain of P, and this proves our claim in this case. 
Note that we can put any power of P there and get the same conclusion. 

Assume now that u solves the wave equation near U only, as in the theorem. Let be as 
above. Fix G U and let V := {(t,x); \t — t%\ < e, \x — x%\ < e} with some < e <C 1. 
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By the non-sharp finite speed of propagation result in [32], u\y depends on u(to, •), to := t\ — 2e, 
restricted to a neighborhood of x\ of size 0(e), assuming that all characteristic cones for the sound 
speed Co := maxc with vertices in V and above t = to, lie entirely in U. The latter is true of 
e< 1. Therefore, we can cut u(to, •) appropriately to make the support compact and contained 
in U fl {t = to}, and then use it as initial data at t = t± — 2e. Then we write v(i, •) = e i (*~*°-' p g, 
where g is the localized Cauchy data at T = to as above. If the cutoff is chosen appropriately, then 
u = v on V and we can apply the argument of the preceding paragraph. 

We continue with the observation that if the set {|t — to I + dist(xo,x) < To} for some to, To, does 
not intersect R x Qq, then one can apply the "smooth version" of the theorem to conclude that 
any solution in a neighborhood of that set that vanishes near [to — To, to + To] x {xo} a l so vanishes 
in that set. 

Next, note that it is enough to show that u = at one of the "extreme" points of U. More 
precisely, let y be such that dist(xo,y) = T. Assume that we can prove that u = near (t,x) = 
(0,y). Let (t, x) be any point in U, and consider the set 

(6.3) |t — t| + dist(xo, x) < T := dist(xo, x). 

It is included in U, and if we know to prove the result for the "extreme" points, then we would get 
that u = near (t, x). 

We divide the rest of the proof into several steps. Let Qq be interior of the complement of the 
maximal unbounded connected component of R" \ T. Let To = d£lo C T. Given x, y both in 
R n \ Qq, let disto(x,y) be the infimum of the length of all smooth curves in R n \ Qq, connecting 
x and y, intersecting T transversally at each common point, where the length is measured in the 
metric c~ 2 dx 2 . 

Choose and fix y E To so that 

(6.4) £:=disto(y,x ) < T. 

If there is no such y, then the theorem follows by its "smooth version". Note that dist(y, xq) < 
disto(y,xo), therefore 

(6.5) [-T + £, T - i] x {y} d U. 

Fix < <5< 1. Let ys be the point in R n \ $7o that lies on the normal geodesic through y, and is 
at distance 5 from y. Let [0, 1] 9 1 1— > j(t) be a smooth curve in R n \ £Iq connecting xq = 7(0) and 
ys = 7(1) of length (in the metric c~ 2 dx 2 ) is n °t exceeding I + 35. It is easy to see that this can 



be done. Assume also that 5 is so small that Is < T, see (6.4). Let < e = disto(7, T)/2. 

Let xi, X2, • • • Xfc be finitely many points on 7, corresponding to increasing values of t, so that 
each segment of c(t) with endpoints (xo,xi), (xi,X2), . . • , (xk,ys) is of length, that we denote by 
V-xqxii e tc, at most e. Apply the "smooth" version of the theorem for the set 

U t, Xl '■= K^); +dist (x,x ) <4 xi} 

with t such that the interval {t; |t — 1| < i Xo xi} is included in [-T, T]. Such t clearly exist and are 
described by |t| < T—£ X0X1 . Next, for any t as above, we have iTj C C/n(RxR n \j7). Therefore, we 
can apply the "smooth" version of the theorem to conclude that u = in a neighborhood of (t, xi) 
for any such t. In other words, u = in a neighborhood of [— T + £ Xo xijT — £ Xo xi] x {xi}- We can 
repeat those arguments to get u = in a neighborhood of [—T+£ XoXl +£ XlX2 ,T—£ XlX2 —£ XlX2 ] x{x2}, 
etc. Thus we get 

n = near [-T + £ s , T - £5} X {y s }. 
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Take the limit 5 — > now to conclude that for any £' > £, u = on an open set containing 
[— T + £', T — £'] x {x; < disto(x, y) <C 1}. Perturb y, and use the fact that our assumptions on 
T are open, i.e., we can perturb T a bit as well, to get 

(6.6) u = in an one-sided neighborhood of [— T + £, T — £] x {y} in R x R™ \ 

and is a solution there. By the assumed regularity of u, taking the trace of u and its normal 
derivative on T, from outside (i.e., in R n \ Qq) is we h defined; and this trace is zero near [— T + 
£, T—l] x {y}. By the transmission conditions ( |3.5[ ), the interior traces of u and its normal derivative 
vanish there as well. Let now c be any smooth extension of c\n to R n . Then 

(6.7) (<9 4 2 - c 2 (x)A)u = in a two-sided neighborhood of [-T + £, T - £} x {y} in R x R". 

By a two-sided neighborhood, we mean a normal one (an open set) — we only used that term to 
emphasize the difference with the neighborhood in ( |6.6| ). Since u vanishes in a neighborhood of 
[— T + £, T — £] x {y} in the exterior, we apply the "smooth" version of the theorem for the speed c 
near sets of the kind [t — 8, t + S] x {y}, < <5 <C 1, \t\ < T — £ — 5 to conclude that u also vanishes 
in the neighborhood in (6.7), after we shrink it if needed by 0(S). We can use our freedom to vary 
T a bit to conclude that 

(6.8) u = in a two-sided neighborhood of [-T + £,T — £] x {y} in R x R n . 



If small enough, that neighborhood is in U, see (6.5). Notice that we actually proved this property 
for any y G R n \ Qq with £ as in (6.4), not only for y G T but the latter case only requires the use 
of the transmission condition. 

For a fixed a 6 Slo, close enough to Tq, let Q± be the maximal open connected component of 
R n \ r containing a. In Ox, we define a distance function disti by minimizing over all smooth 
curves, transversal to T, that stay in Let T\ = T — £ > 0, see (6.4), and choose y\ G T, if 
possible, so that 

(6.9) h := disti (yi,y) <T X . 

Use the arguments above with T replaced by T\, xq replaced by y before, and y replacing y\. The 
only difference is that we work in £l\ now. We then get 

(6.10) u = in a two-sided neighborhood of [-T + £ + £ 1: T - £ - l x ] x {yi} in R x R". 

Let now z be a point with dist(z,a?o) < T. Choose 5 > so that there is a smooth curve 7 of 
length dist(z, Xq) + 5 crossing V transversally each time, and connecting xq and z. Let also <5 ^ 1 
so that dist(z, xq) + 5 <T. By a compactness argument, 7 will cross T finitely many times. Apply 
the argument above for each segment either in Oo or in R" \ £Iq to get 

u = near [— T — dist(z, xq) — S,T + dist(z, xq) + 5] x {z}. 

This shows that u = in the interior of U. Since our assumptions allow us to increase T slightly, 
u = near U as claimed. □ 



7. WE CAN GET THE WHOLE CAUCHY DATA ON [0, T] X <9$7 

We will show here that knowing Af , one can recover the Neumann derivative of the solution at 
[0, T] x dQ as well. This is done by applying a non-local ^DO to Af . This is known in principle, 
and we do not use it in our proofs. It reveals links to Control Theory however. 
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We will define first the outgoing DN map. Given g G Cq°([0, oo) x an), let w solve the exterior 
mixed problem with c = 1: 

in(0,T)xR n , 

9, 
0. 
0. 



(7.1) 



t ~ A)w 
w \[o,T]xdn 
w\ t =o 
dtw\ t =o 



Then we set 



Ng 



dw 
dv 



[0,T]xdn 



By [T7], for g G iTL([0,T] x 90), we have [w,w t ] G C([0,T); H); therefore, 

n : ^([cr] x an) -> c([o,t] x #3 (an)) 

is continuous. Note that the results in [T7] require the domain to be bounded but by finite domain 
of dependence we can remove that restriction in our case. We also refer to [H Proposition 2] for a 
sharp domain of dependence result for exterior problems. 



Lemma 7.1. Let u solve (2.1) with f G % compactly supported in n. Assume that c = 1 outside 



n. T/ien /or any T > 0, Af determines uniquely u in [0, T] x R n \ n and t/ie normal derivative of 
u on [0, T] x <9n as follows: 

(a) The solution u in [0, T] x R n \ n coincides with the solution of (7.1) with g = Af, 

(b) We have 

NAf. 



(7.2) 



dv 



[0,T]x<9fi 

Proof. Let u> be the solution of (7.1) with g = Af. The latter is in Hh-,([Q, T] x an), see the 



paragraph after (2.2). Let u be the solution of (2.1). Then u — w solves the unit speed wave 
equation in [0, T] x R n \ n with zero Dirichlet data and zero initial data. Therefore, u = w in 
[0, T] x R n \ n. □ 

Remark 7.1. Note that c = 1 outside n was not a necessary assumption. 
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